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A  COHERENT  NONLINEAR  THEORY  OF  AURORAL 
KILOMETRIC  RADIATION 


I .  INTRODUCTION 

There  have  been  several  observational  studies  made  in  recent  years 
of  high  intensity  kilometric  radiation  of  terrestrial  origin.  Measure¬ 
ments  have  been  made  by  satellite  from  the  source  region,  in  the  auroral 
zone  at  R  2-3  R^,  all  of  the  way  out  to  R  25-30 
the  source  region. 

A  number  of  interesting  features  of  the  auroral  kilometric  radia¬ 
tion  (AKR)  have  been  deduced  as  a  result  of  these  satellite  measure¬ 
ments.  The  first  measurements  were  made  on  Ogo  I  [Dunkei,  et  al . , 

1970].  More  detailed  measurements  were  made  on  Imp  6  and  Imp  8,  which 
were  analyzed  by  Garnett,  et  al.  [1974].  Gurnett  found  the  spectrum  to 
lie  in  50  kHz  <  f  <  500  kHz,  with  peak  intensity  at  f  200  kHz.  The 
peak  intensity  power  emission  was  estimated  to  be  about  10^  W.  He  noted 
that  the  radiation  appears  to  originate  at  low  altitudes  (R  2-3  R^) 
in  the  auroral  region,  and  to  be  closely  associated  with  the  occurrence 
of  discrete  auroral  arcs,  which  are  believed  to  be  generated  by  intense 
inverted  V  electron  precipitation  bands.  When  arcs  do  not  occur,  only 
a  small  band  of  diffuse  aurora  is  present,  and  the  radiation  disappears. 
The  radiated  power  was  estimated  to  be  close  to  1  percent  of  the  maximum 
energy  dissipated  by  the  auroral  charged  particle  precipitation. 

Recent  observations  in  the  source  region  from  Isis  I  reveal  and 

establish  further  properties  of  the  AKR,  which  were  not  well-established 

from  the  far-field  measurements  [Benson  and  Calvert,  1979] .  Those 
Note:  Manuscript  aubmitted  August  1,  1979. 
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far  away  from 
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observations  showed  that  the  radiation  was  generated  in  the  X-mode 
just  above  the  local  cutoff  frequency,  and  propagated  almost  perpen¬ 
dicular  to  the  background  magnetic  field.  It  was  found  to  be  generated 

within  density-depleted  regions  with  peak  density  such  that  a)  <0.2 

pe 

u)  .  (These  are  probably  caused  by  the  large  potential  drop  in  those 
ce 

regions  associated  with  the  dis-rete  auroral  arcs.) 

A  number  of  theories  have  been  proposed  in  the  literature  to 
explain  the  occurrence  of  the  high-intensity  radiation.  Oya  (19741 
proposed  a  mechanism  for  the  well-known  Jovian  decametric  radiation 
(thought  to  have  a  similar  origin  as  that  of  the  terrestrial  kilometric 
radiation)  which  involved  the  creation  of  electromagnetic  radiation  by 
a  mode  conversion  process.  This  mechanism  was  extended  to  kilometric 
radiation  by  Benson  [19751.  In  this  model,  excited  plasma  waves  were 
converted  into  the  X-mode  at  the  upper  hybrid  frequency.  This  mode  was 
then  converted  to  the  0-mode.  The  overall  efficiency  of  this  mechanism 
is  estimated  to  be  a  few  percent. 

There  are  two  major  drawbacks  with  this  theory.  First,  the  polar¬ 
ization  ia  in  disagreement  with  recent  polarization  measurements, 
including  both  far-field  measurements  [Gurnett  and  Green,  19781  and 
measurements  at  the  source  [Benson  and  Calvert,  19791 .  Second,  the 
overall  efficiency,  which  is  the  product  of  the  efficiencies  of  two 
inefficient  conversion  mechanisms  is  too  small.  This  does  not  even 
take  into  account  the  conversion  efficiency  of  precipitating  electrons 
into  electrostatic  plasma  waves  in  the  first  place.  Thus  it  is  highly 
unlikely  that  this  indirect  mechanism  could  produce  final  output 


radiation  which  is  on  the  order  of  1  percent  of  precipitating  electron 
(or  electron  beam)  energy. 

Melrose  [1976]  proposed  a  theory  for  the  Jovian  decametric  and 
terrestrial  kilometric  radiation  which  involves  directly  amplified  gyro- 
emissions.  In  this  mechanism,  radiation  is  emitted  at  the  electron 
cyclotron  frequency  and  its  harmonics,  which  have  been  Doppler  Shifted 
by  electron  streaming.  This  radiation  is  then  amplified  due  to  assumed 
anisotropic  velocity  distribution  functions.  It  is  also  assumed  that 
the  electron  streaming  velocity  is  large  enough  to  allow  the  radiation  to 
escape  to  free  space. 

The  frequency  for  the  radiation  required  in  Melrose's  theory  will 
turn  out  to  be  the  most  favorable  frequency  for  radiation  in  the  theory 
proposed  in  the  present  paper.  However,  to  get  growth  rates  required  to 
produce  kilometric  radiation  at  the  observed  level  by  anisotropic  velocity 

space  instabilities  would  appear  to  require  very  large  anisotropies  in 
velocity  space.  There  is  at  present  no  experimental  evidence  for  such 
large  anisotropies. 

Wu  and  Lee  [1979]  have  recently  proposed  a  mechanism  which  involves 
a  velocity-space  instability,  namely  a  loss-cone  type  instability.  In 
their  model,  precipitating  electrons  stream  down  the  (converging) 
magnetic  field  lines  toward  Earth  until  they  encounter  an  effective 
magnetic  mirror  in  a  region  of  higher  magnetic  field  strength.  They 
then  undergo  reflection  and  travel  up  diverging  magnetic  field  lines. 

Some  of  the  upgoing  particles  are  lost  through  the  loss  cone,  and  the 
resulting  velocity  distributions  lead  to  a  loss  cone  instability. 
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Wu  and  Lee  incorporate  into  their  model  a  long  density-depleted 
cavity  in  which  the  unstable  wave  grows,  as  found  by  Benson  and  Calvert 
on  Isis  I.  This  cavity  was  not  considered  in  any  of  the  previous 
papers,  but  will  be  assumed  in  our  theory.  The  most  probable  explana¬ 
tion  is  that  the  cavity  is  formed  by  large  potential  drops  known  to 
be  present  in  the  auroral  region,  which  also  tend  to  drive  large 
electron  (and  even  ion)  beams.  For  example,  if  we  assume  the  back¬ 
ground  density  is  close  to  thermodynamic  equilibrium,  then  the  density 
profile  is  determined  by  the  Boltzmann  form 

n  (x)  =  n  exp  [-  e<ti(x)/KT  ] 
e  o  e 

where  v{)(x)  is  the  local  potential  and  T^  the  electron  temperature. 

Thus  the  regions  of  high  potential  have  lower  density. 

One  drawback  of  the  Wu-Lee  theory  is  that  to  get  the  large  growth 

lengths  required,  loss-cone  angles  of  9  >  45-60®  appear  to  be  required, 

as  for  example,  the  effective  angles  in  the  distribution  functions 

assumed  by  Wu  and  Lee.  The  observed  loss  cone  angles  in  the  region  of 

origin  are  more  like  9  '''  15-20®. 

Another  recent  theory  for  kilometric  radiation  was  proposed  by 

Roux  and  Pellat  (1979) .  Their  mechanism  involved  the  nonlinear 

beating  of  electrostatic  waves  near  the  upper  hybrid  to  produce  an 

electromagnetic  wave.  Their  theory  predicts  that  the  0-mode  should 

dominate  at  ui  =  w  ^  (at  the  upper  hybrid  frequency)  for  u)  >  w  , 
uh  j  ^  j  pg 

and  that  the  X-mode  should  dominate  at  uJ  =  2  a)  ,  for  OJ  <  w  . 

uh  pe  ce 


Roux  and  Pellat  argue  that  a  coherent  nonlinear  theory,  such  as 
their  proposal  or  the  present  proposal,  should  be  more  efficient  in 
general  than  an  incoherent  linear  mechanism  or  even  a  coherent  linear 
mechanism.  However,  their  particular  nonlinear  mechanism,  which  re¬ 
quires  the  creation  of  a  high  phase  velocity  transverse  wave  from  two 
low  phase  velocity  longitudinal  waves,  would  be  expected  to  be  rather 
inefficient.  Also,  their  prediction  of  radiation  in  the  X-mode  at 
w  =  the  density  depleted  case  observed  is  in  apparent  dis¬ 

agreement  with  the  observation  that  the  frequency  at  the  origin  is  just 
above  the  right-hand  cutoff  frequency.  Roux  and  Pellat  also  mention 
the  possibility  of  the  nonlinear  interaction  of  an  upper  hybrid  wave 
and  a  lower  hybrid  wave  to  produce  radiation  at  U)  =  However, 

since  w  <<  to  in  the  density  cavity  then  this  would  give  to  =  1.02  to  , 
pe  ce  ce 

and  the  resulting  radiation  may  be  below  the  right  hand  cutoff.  Also, 
the  efficiencies  for  this  case  would  appear  to  be  less  than  the  case 
for  two  upper  hybrid  waves,  but  it  has  not  yet  been  analyzed. 

Finally,  a  theory  proposed  by  Palmadesso ,  et  al.  [1975],  assumed 
that  electromagnetic  "noise"  gets  amplified  through  its  interaction 
with  ion  turbulence.  The  nonlinear  beat  wave  of  these  two  waves  could 
be  amplified  by  beams,  and  thus  the  electromagnetic  wave  would  be 
paraunetrically  amplified.  We  will  employ  a  similar  mechanism  in  the 
theory  presented  in  the  present  paper. 

In  the  Palmadesso  theory,  both  the  RX  and  LO  modes  were  assumed  to 
be  amplified  by  the  above  mechanism.  However,  for  the  beam  to  signifi¬ 
cantly  amplify  the  beat  wave,  and  hence  the  electromagnetic  wave,  the 


beat  wave  was  assumed  to  be  almost  a  natural  mode  of  the  plasma,  and 
required  to  have  a  slow  phase  velocity  (w/k  <<  c) .  This  required  gene¬ 
ration  at  frequencies  below  the  upper  hybrid.  Thus  the  RX  mode  would 
encounter  a  cutoff  at  frequencies  between  the  upper  hybrid  and  the 
right-hand  cutoff  and  not  escape  to  free  space.  Thus,  only  the  LO 
mode  would  be  observed,  contrary  to  recent  polarization  measurements. 
Also,  the  ion  waves  were  taken  to  be  incoherent  (turbulent) ,  and 
assumed  to  be  isotropic  for  convenience  of  calculation.  Recent  measure¬ 
ments  of  EIC  waves,  e.g.,  Temerin,  et  al .  [1978],  show  strong  electro¬ 
static  ion  cyclotron  (EIC)  wave  fluctuations  which  are  coherent  and 
primarily  perpendicular  to  the  background  magnetic  field.  The  coherence 
of  the  waves  was  unanticipated  and  suggests  a  modified  application  of 
the  mechanism  in  Palmadesso,  et  al . ,  can  be  employed. 

We  will  construct  the  present  theory  with  the  same  principle 
underlying  mechanism.  However,  we  will  construct  a  new  theoretical 
framework  so  the  formulation  of  the  theory  will  be  different.  We  will 
assume  coherent  electrostatic  density  fluctuations  (due  to  low- 
frequency  EIC  waves)  perpendicular  to  the  magnetic  field.  We  will 
concentrate  on  the  case  that  there  is  a  resonant  interaction  between 
electromagnetic  wave  and  the  coherent  density  fluctuations,  a  case  not 
considered  by  Palmadesso ,  et  aJ .  It  will  be  shown  that  the  X-mode 
is  preferentially  amplified  over  the  0-mode  in  this  model,  and  that 
the  X  mode  can  grow  rapidly  in  a  narrow  frequency  band  that  is 
accessible  to  free  space.  The  basic  requirement  for  this  last  condi¬ 
tion  to  exist  will  be  shown  to  be  that  a  density  cavity  of  precisely 
the  density  observed  is  present  at  the  source. 
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In  Sec.  II  we  will  introduce  the  basic  model  for  the  AKR,  derive 


the  wave  equation  and  associated  dispersion  relation,  and  discuss  solu¬ 
tions.  In  Sec.  Ill  we  will  analyze  the  wave  equation  to  find  the  in¬ 
stability  criterion  and  the  growth  rate  for  the  unstable  modes.  In 
Sec.  IV  we  discuss  in  more  detail  the  physics  of  the  wave  amplification. 


Results  are  summarized  in  Sec.  V. 


II.  MODEL  AND  GEOMETRY 


We  assume  that  initially  we  have  electromagnetic  radiation  propa¬ 
gating  at  the  noise  level  in  the  auroral  zone,  which  has  been  generated 
by  one  of  a  variety  of  possible  mechanisms:  Cerenkov,  electron  cyclo¬ 
tron,  bremstrahlung ,  etc.  We  also  assume  that  there  are  coherent  density 
fluctuations  perpendicular  to  the  B-field  due  to  EIC  waves,  as  observed 
by  Temerin  and  coworkers.  The  background  magnetic  field  =  B^z  is  taken 
along  the  z  direction.  The  electromagnetic  wave  is  assumed  to  propagate 
almost  perpendicular  to  the  magnetic  field,  with  propagation  vector 
k  =  (k  ,  0,  k  ),  where  k-  <<  k^ ,  as  observed  by  Benson  and  Calvert 

.X  Z  Z  X 

(19791 .  (See  Fig.  1.) 

There  are  two  independent  modes  of  the  electromagnetic  wave,  with 
E-field  vectors: 


E  = 

[ExIx)^ 

>  Ey(x)?] 

likyZ  -  Lijt) 
e  ^ 

X-mode 

(la) 

E  = 

|^Ex(x)x 

+  E^(x)zJ 

i  (k,2  -  uJt) 

e  * 

0-mode 

(lb) 

We  are  dealing  with  high-frequency  waves  this  model  so  that  cj  ~  u)  , 

ce 

uj  .  The  waves  are  primarily  transverse  waves,  so  E  <<  E  ,  E  . 
pe  ^  X  y  z 

The  EIC  waves  have  associated  density  fluctuations  of  both  ions 
and  electrons:  the  electrons  are  assumed  to  follow  the  ions  in  the 
fluctuations  so  that  the  resulting  total  plasma  fluctuation  is  quasi¬ 
neutral.  Then  the  plasma  density  is  of  the  form 

n(x,t)  =  n  +  n, (x,t)  (2a) 

o  * 

nj(x,t)  =  ,6n  cos  (k^x  -  uj^t)  (2b) 


8 


where  is  the  background  density,  nj  is  the  fluctuating  density,  and 

u)^  and  k^  are  the  frequency  and  the  wavenumber  of  the  EIC  wave.  Since 

u).  A-  u)  .  <<  u)  Xu),  the  electromagnetic  wave  sees  the  density  fluctu- 
1  cl  ce 

ation  as  stationary.  Thus  we  can  ignore  the  oj^t  term  in  nj  in  our 
electromagnetic  wave  equation: 


CO 


(3) 


where  and  are  the  dielectric  constants  perpendicular  to  and 

along  the  magnetic  field,  respectively,  and  K  is  the  Hall  term. 

H 

We  will  divide  the  components  of  K  according  to 


K.  .  =  K.  .  +  <5k.  .  (5) 

1]  1^0  ID 


where  K  .  are  the  tensor  components  for  the  case  of  n(x,t)  =  n  ,  and 
130  o 

5k, ,  are  the  perturbations  introduced  in  by  the  density  flvrrtuation 


n^^.  Including  the  effect  of  a  beam  of  velocity  =  v^z  along  the 
B-field,  the  dielectric  tensor  components  for  the  nonf luctuatinq 


f 


background  plasma  for  the  case  of  a  cold  plasma  and  a  warm  beam  are: 


10  ^  (0^ 

V  =  1  _  _  _ b 

J.O  co^-to  ^  ((0-k  V  +iv)^-(o  ^ 

ce  zb  ce 


10 


ce  I  peo 

K„_  =  —  + 


Ho  (0 


(0 

C0^-(0 


(Of 


((0-k  v^+iv)^-(o  ^ 
ce  zb  ce 


2  2 

(0  (O^ 

K  =  1  .  -E®° _ ^ _ - 

Ho  (0"^  ((0-k  v,+iv)^ 

z  b 


(6a) 


(6b) 


(6c) 


with  the  plasma  frequency  for  plasma  density  n^.  We  have  neglected 

the  low-frequency  ion  terms,  since  they  are  not  of  interest  for  the 
high-frequency  electromagnetic  wave.  Here  (o^^  is  the  beam-plasma  fre¬ 
quency,  and  V  enters  because  of  the  finite  beam  width  themal  effects 
for  an  assumed  Lorentzian  beam  velocity  profile.  In  terms  of  the  beam 
spread  Av  in  velocity  space, 

V  =  k  Av  (7) 

z 

The  use  of  a  Lorentzian  beam  profile  for  thermal  effects  is  valid 
provided  v^^  »  v^  Av  where  v^  is  the  average  electron  thermal  velocity 
of  the  beam. 

From  the  form  of  E  for  the  X  and  0  mode  given  by  Eg.  (1) 
we  may  obtain  a  reduced  wave  equation  for  each  mode  through  Eq.  (3) . 

92e 


“  ^Ci2  k^x)  E^  =  0 


z 

1^ 


(6i  +  eS2  cos  k^x)  E^  =  0 


X-mode 


0-mode 


(8a) 


(8b) 
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where 


/  \^  K  ^  ^ 

~  /  0)  \  ^  Ho  2 

\=/  [_  J 

UJ^U)  2  r  2  kJ 

^  peo  ,  HO 

cMur-aT^)  (K,  -n^)" 

ce  lo  z 


(K,  -n-^)K,, 
lo  Z  Ho 


^  ^  (k7-n^  ^ 

lo  z 


^Ho 

(K,  -n  -) 
lo  z 


/w^  \  cj^/K  -n^ 

"  ^.o  (^)  ■>;  -  ^ 


The  form  of  Eq.  (8)  is  obtained  under  the  assumptions; 


1 10^  -  CO  ?  I  >''  £  CO  ^ 


e-  <<  1 


(10a) 


Since  measured  density  fluctuations  give  £  <  0.5,  the  £“  <<  1  assumption 
appears  to  be  quite  good.  Also,  Eq.  (10a)  is  quite  valid  except  for  a 
very  narrow  range  around  the  upper  hybrid  frequency: 


1.01  (0  ^  u)  0.99  10 

uh  uh 


The  latter  frequency  range  will  not  be  of  interest  to  us  because  it 
lies  below  the  right  hand  cutoff,  and  involves  a  wave  inaccessible  to 
free  space. 

Thus  the  wave  equations  for  both  the  X  and  the  O  mode  reduce  to 


Mathieu's  equation  for  the  model  being  used  here.  From  this  equation 


i 


we  may  find  stability  conditions  and  growth  rates  which  arise  from 
resonant  spatial  interaction  between  the  electromagnetic  wave  and 
spatial  density  fluctuations,  as  well  as  the  modification  of  the 
electromagnetic  waves  by  the  density  fluctuations. 

For  purposes  of  illustration  of  the  solution,  we  will  solve  the 
wave  equation  for  the  X-mode.  The  0-mode  solution  is  quite  similar. 
In  obtaining  Eq.  (8a),  we  used  the  relationship  between  and  E^ 
from  Eq.  (3).  From  that,  the  appropriate  form  of  the  solution  for  E 
of  the  X-mode  is 


E  =  A'  y  - 


Ho 


i  (k_z-u)t) 


(K,  -n  -) 
lo  z 


±  ±i(q+nk. 


.  )X 


(12) 


where  a  ~  =  1,  and  for  no  density  fluctuations 
o 


a  (e  =  0)  =  0  (n  0) 
n 


(13) 


By  substituting  this  form  into  our  wave  equation,  Eq.  (3) ,  we  find 


the  dispersion  relation  to  order 


I  2  2 


qS  -  q'^  (2k?  +  3ai)  -  q^  I -y*"  - 


4  2\ 

i  -  3ai  j 


[2ai (ai-k^ 


)2  +  (k?  -  ai)J  +  0{c‘*)  =  0 


(14) 


Also,  solution  of  the  coefficients  to  order  e  yields 

*  -ea. 


“  2Iai+(q-k^)2] 


(15a) 


(15b) 


f 


t  2 

"  2[ai+(q-k.)^] 

These  equations  show  the  interaction  between  the  electromagnetic  wave 

and  the  first  harmonic  of  the  density  fluctuation.  If  we  included  in 

higher  orders  of  e,  we  would  get  the  higher  harmonics.  In  general, 

for  the  inclusion  of  harmonics  up  to  n,  we  get  a  dispersion  relation 

of  order  2n+l  in  q^,  and  coefficients  a  ~  .  It  should  be  noted 

n 

that  in  the  limit  e  -►  0,  the  dispersion  relation  in  Eq.  (14)  has  roots 


‘31,2 

= 

(16a) 

<33,4 

= 

-  k^) 

(16b) 

q5,6 

= 

(t/oT  +  l^i) 

(16c) 

The  first 

two 

roots  are  just 

the  fundaunental  mode,  and  the  last 

four 

are  simple 

first 

harmonics. 

However,  for  £  -►  0,  the  amplitudes 

+ 

of  these  harmonics  are  zero. 
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III.  STABILITY  CRITERION  AND  GROWTH  RATES 


When  the  electromagnetic  wave  is  propagating  (i.e.,  for  aj  '‘0), 
then  the  q's  in  Eq.  (16)  for  the  unmodulated  wave  are  real.  However, 
the  presence  of  density  fluctuations  (so  e  ^  0)  causes  the  roots  for  q 
in  Eq.  (14)  to  be  complex  for  certain  conditions,  so  that  the  propaga¬ 
tion  wave  has  growing  or  decaying  solutions.  The  general  requirement 
for  such  an  instability  to  exist  is 

»/a7  =  pk^  (X-mode),  (17a) 

=  pk^  (0-mode)  (17b) 


where  p  is  a  positive  integer.  The  larger  6  is,  the  less  strictly 
this  criterion  has  to  be  satisfied.  A  stability  diagram  is  shown  in 
Fig.  2.  Basically  Eq.  (17)  is  just  the  condition  for  resonant  inter¬ 
action  between  the  electromagnetic  wave  and  a  spatial  harmonic  of 
the  density  fluctuation.  As  one  would  expect,  the  p  -  1  model  is  the 
fastest  growing. 

The  conditions  for  instability  of  the  X-mode  may  be  rewritten 


.2 


^  (K  -n^) 
1  z 


(18) 


Generally,  ck^/u)  ~  l  since  the  ion  wave  wavelength  tends  to  be 
shorter  than  the  free-space  electromagnetic  wavelength  (but  not  of 
the  electromagnetic  wavelength  in  the  pldsma  medium) .  Thus  for  there 
to  be  a  solution  of  the  X-mode  stability  criterion,  we  must  have 


'1 

OJ  “ 

ce 


(co-k  V  )  ■ 
z  b 


u)  *■ 
pe 

>  ^  ' 


ce 


(19) 


\ 
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There  are  two  frequency  regimes  where  this  is  satisfied; 

<i)  u)  ~  ^ce"  frequency  regime  has  waves  which  are  not 

accessible  to  free  space,  because  of  the  cutoff  layer  between  the 
upper  hybrid  and  right-hand  cutoff.  Also,  there  may  be  electron 
cyclotron  damping  which  would  partially  offset  any  growth  of  the  wave. 
Thus  this  regime  is  unlikely  to  produce  high- intensity  radiation  which 
is  observable  in  free  space. 

(ii)  u)  *<  Doppler-shifted  beam-cyclotron 

resonance.  There  are  waves  in  this  regime  which  are  accessible  to  free 
space  provided 


CO  /  CO 

U)  ^  k  V^  >  <0  =  >  CO  2  + 

ce  zb  R  2  \  pe  4 


Expanding  the  right  hand  side  for  co  2  <<  qj  2  this  becomes 

pe  ce 

CO  2  <  k  V,  (0 
pe  z  b  ce 


For  excimple,  for  k  0.3k  and  v  c/4,  we  find 
z  b 


(0  <  0.  3  CO 

pe 


This  is  to  be  compared  to  co  <  0.2  co  in  the  source  of  AKR.  Thus  the 

pe 

density  depletion  observed  by  Isis  I  at  the  source  is  of  just 

the  amount  predicted  by  theory  to  be  necessary  for  radiation.  When 

the  condition  in  Eq.  (22)  is  satisfied,  we  get  radiation  in  the  narrow 


CO  +  CO  2/(jj  <  (0  <  CO  +  k  V, 

ce  pe  ce  ce  z  b 


T 


just  above  the  right  hand  cutoff.  This  frequency  regime  is  then  the 
most  probable  regime  of  the  kilometric  radiation.  It  should  be  noted 
that  there  ’s  a  range  of  values  of  the  plasma  and  cyclotron  frequencies 
at  the  source  of  AKR  in  the  aurora,  and  the  frequencies  which  lie  in 
the  narrow  band  given  by  Eq.  (23)  are  different  for  each  sourcepoint. 
Thus  the  total  radiation  received  has  the  broad  spectrum  between  50 
and  500  kHz  observed.  Also,  Eq.  (21)  shows  that  radiation  will  not 
be  observed  in  regions  where  the  beeuns  are  not  of  sufficiently  high 
energy,  or  where  the  plasma  density  is  too  high  (such  as  in  regions 
R  <  2  R^,  or  outside  of  the  density  depleted  region) ,  or  where  the 
magnetic  field  is  too  low  (such  as  in  regions  at  R  >  3  R^) . 

Having  established  the  frequency  regime  for  the  radiation,  it 

must  also  be  pointed  out  Eq.  (19)  also  gives  a  minimum  required  beam 

density  n  ,  namely 
b 


i)  ■  ii)  ■ 


(24) 


where  we  have  included  in  the  finite  beam  width  and  damping  term  v. 
On  resonance,  this  reduces  to 


(25) 


which,  for  \)  =  kAv'v^0.3kv,  k  "^O.lk,  v  ^O.lc  gives 
2  z  b  z  b 

n,  '*  1.5  ><  10‘^  n 
L  o 


("le) 
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Recall  is  the  reduced  plasma  density  due  to  density  depletion,  so  in 
this  case  instability  would  be  possible  for  ^  ^  ^  10"^  n^,  where  n^ 
is  the  background  plasma  density  outside  the  density  cavity. 

For  the  0-mode  we  may  make  the  same  analysis  as  we  did  for  the 
X-mode.  For  instability  for  this  mode  we  require 


Now  the  left  hand  side  is  always  greater  than  one,  whereas  the  right 

hand  side  is  always  less  than  one  for  propagation  waves.  Thus  we  have 

no  unstable  solution  of  this  equation  and  high-intensity  radiation  is 

not  possible.  It  should  be  recalled  that  Eg.  (27)  breaks  down  in  a 

very  narrow  band  near  the  upper  hybrid.  Thus  the  presence  of  a  k^  /  0 

makes  it  conceivable  that  there  might  be  an  instability  and  some 

radiation  in  the  0-mode  near  co  =  o)  ,  .  However,  even  if  such  an  in- 

uh 

stability  can  exist,  the  band  is  very  narrow,  and  the  wave  equation  is 
not  nearly  as  unstable  as  the  Mathieu  equation.  Thus,  in  any  case,  the 
X-mode  radiation  would  definitely  dominate,  and  we  will  concentrate  on 
that  mode  from  now  on. 

We  may  now  calculate  the  growth  rate  for  our  X-mode  radiation. 

For  excunple,  the  dispersion  relation  for  the  wave  equation  including 
the  first  harmonic  interactions,  Eq.  (14)  ,  will  yield  an  instability 
provided  the  equation  has  a  complex  (non-real)  solution.  The  spatial 
growth  rate  for  the  first  harmonic  instability  which  occurs  when 
/a7  «  k^  (p  =  1  mode)  is  for  e  «  1 


I 
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be  1000-2000  tan,  enough  to  allow  the  wave  to  propagate  to  a  region 
where  the  background  density  is  such  that  at  the  cavity  boundary  is 
then  accessible  to  free  space.  Note  that  wave  propagating  in  directions 
toward  higher  magnetic  fields  and  higher  background  density  will 
eventually  reach  the  right  hand  cutoff  layer  and  reflect  back  toward 
lower  field  and  density  regions.  Thus  this  radiation  may  also  reach  a 
region  of  escape  to  free  space.  In  any  case,  the  path  length  avail¬ 
able  is  quite  adequate  to  produce  the  growth  required. 
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IV.  PHYSICAL  MECHANISM 

Now  that  we  have  found  the  unstable  solutions  to  the  wave  equation 
and  shown  that  they  are  adequate  to  explain  the  kilometric  radiation 
observed,  it  is  informative  to  look  at  the  physics  of  the  amplification 
in  more  detail . 

As  shown  in  Section  III,  there  were  two  frequency  regimes  in  which 

growing  waves  could  exist;  u)  ~  u)  and  u)  ~  u)  +  k  v.  .  The  former  was 

ce  ce  z  b 

discounted  because  of  inaccessibility  and  possible  cyclotron  damping 
problems.  Let  us  consider  the  energy  density  of  the  wave  in  the  two 
frequency  ranges  as  seen  by  the  electrons  in  their  rotating  frame: 


(30) 


where 


K' 

o 


0)  “ 

1  - 


(ui '  -k  V,  ) ' 

z  b 


(31) 


and  a)'  =  u)  -  *  Doppler  shifted  frequency  seen  by  rotating  electron. 


£ 

u.f 


r/  ' 

(  1  + 

(u'+k  V  )\ 

b  zb 

ip  I  ^ 

Ip  |2  \ 

l\  “ 

(w'-k  V  ) 3  j 
zb  / 

) 

(32) 


Now  U  is  negative  when  u)'  <  k  v,  ,  i.e.,  near  the  Doppler  shifted  cyclotron 

z  b 

resonance  a)  ~  u)  +  k  v,  the  X-mode  is  a  negative  energy  mode.  However, 
ce  z  b 

for  oj  <“  0)^^,  it  is  very  difficult  to  get  a  negative  energy  mode.  Thus 

the  former  mode  can  absorb  energy  from  the  beam,  whereas  the  latter 

cannot,  because  of  the  relation  of  the  phase  velocity  in  the  electron 

frame  u)'/k  to  the  beam  velocity  v,  . 

z  b 


\ 
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When  the  frequency  range  to  <”  is  not  a  negative  energy  mode, 
any  energy  it  absorbs  must  come  from  the  EIC  density  fluctuations  via 


resonant  interaction.  [See  Fig.  8(a).)  Ordinarily,  energy  flow  from 


the  electromagnetic  wave  to  the  lower  energy  ion  waves  would  be  thermo 


dynamically  favored,  and  the  X-mode  is  more  li)cely  to  be  absorbed  in 


energy  flow  from  the  beam  into  the  negative  energy  X-mode  via  coherent 


interaction  with  the  ion  mode  is  thermodynamically  favored.  [See 


In  Fig.  8(b)  it  is  seen  that  the  beam  interacts  with  and  effec 


tively  amplifies  the  nonlinear  beat  wave  between  the  EIC  mode  and  the 


X-mode.  This  mode  is  a  driven  or  virtual  mode 


mode  since  it  is  not  necessarily  a  natural  mode  of  the  system.  We  may 


learn  more  about  this  nonlinearly  driven  mode  by  looking  at  its  phase 


which  is  the  phase  difference  between  the  modes  producing  it.  For 


interaction  with  the  pth  spatial  harmonic 


Thus 


which  is  just  the  phase  of  the  right-hand  circularly  polarized  (R)  mode 


travelling  almost  along  the  magnetic  field.  This  mode  has  an  unstable 


resonance  at  the  Doppler-shifter  beam  cyclotron  resonance  and  is  thus 


r 


amplified  by  the  beam.  The  amplification  of  the  virtual  mode  then 
amplifies  the  parent  waves  which  are  driving  it.  It  should  be  noted 
that  electromagnetic  noise  is  not  even  necessary  initially.  The  beam 
may  drive  an  R-mode,  which  interacts  with  EIC  waves  to  produce  X-mode 
waves.  (See  Fig.  9.) 
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V.  SUMMARY  AND  CONCLUSIONS 

Nonlinear  amplification  of  electromagnetic  wave  through  coherent 
interaction  with  EIC  waves  and  electron  beauns  for  the  model  used  appears 
to  be  quite  adequate  to  produce  the  power  levels  of  kilometric  radiation 
observed.  All  aspects  of  the  model  and  results  appear  to  be  in  good 
agreement  with  current  observations.  The  radiation  is  predicted  to  lie 
in  a  narrow  frequency  band  just  above  the  local  right-hand  cutoff,  as 
observed  at  the  source.  The  basic  conditions  for  this  to  occur  in  our 
theory  is  that  we  have  a  high-energy  beam,  and  a  density-depleted  cavity 
of  precisely  the  type  observed.  For  large  growth  due  to  interaction 
with  density  fluctuations,  the  radiation  must  propagate  almost  perpen¬ 
dicular  to  the  background  magnetic  field  and  have  X-mode  polarization, 
both  characteristics  just  as  observed  by  Isis  I  source  measurements. 

The  existence  of  the  various  p-mode  instabilities  indicates  a  harmonic  band 
structure  consisting  of  the  lowest  p-modes  might  be  observable,  in  agree¬ 
ment  with  recent  measurements  [Shawhan,  1979] . 

For  the  p  =  1  mode,  growth  lengths  of  L  100  km  seem  to  be  adequate 
to  bring  the  kilometric  radiation  out  of  the  noise  to  the  power  levels 
observed.  Imperfections  in  the  source  region  may  raise  this  somewhat. 
Higher  order  modes  may  take  several  hundred  kilometers.  Such  path 
lengths  appear  to  be  easily  obtained. 

The  locally  generated  frequency  should  decrease  as  the  radius  R 

increases.  The  model  implies  that  the  reason  the  radiation  is  restricted 

to  R  't  2-3  R  is  that  for  R  <  2  R  the  plasma  density  is  too  high  and 

£« 

for  R  >  3  the  magnetic  field  is  too  low,  to  put  the  Doppler-shifted 


\ 


boam  cyclotron  frequency  above  cutoff.  Only  for  R  'v  2-3  can  radiation 
bo  qeneratod  above  cutoff. 

It  should  be  noted  that  there  are  possible  future  applications 
of  the  proposed  mechanism  for  generation  of  radiation  to  laboratory 
sources  of  high-intensity  radiation,  and  possibly  even  plasma-based 
masers  or  lasers.  This  is  related  to  the  Kaw,  Lin  and  Dawson's  (1973] 
proposal  for  a  plasma  laser.  A  future  paper  on  that  area  is  planned. 
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unstabi* 


stability  diagram  for  the  wave  equation  of  the  X-mode.  In 
stability  occurs  whenever  -  pk^,  where  p  is  an  integer 
As  the  relative  density  fluctuation  amplitude  e  increases, 
the  less  strictly  this  criterion  must  be  satisfied. 


7.  Path  of  electromagnetic  ray  in  density-depleted  region.  Multiple 
reflections  off  of  cavity  boundaries  may  provide  long  paths 
for  growth.  The  way  may  eventually  reach  a  weak  field  and/or 
low  plasma  density  region  where  it  is  accessible  to  free  space. 
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EM  wave 


EIC  wave 


Fig.  8.  (a)  Simple  resonant  scattering  of  the  EM  and  EIC  waves.  Energy 

transfer  takes  place  in  this  scattering,  normally  from  the  EM 
wave  to  the  EIC  wave.  (b)  Resonant  scattering  of  EM  and  EIC 
waves,  in  which  the  beat  "quasimode"  interacts  with  the  beam 
at  resonance.  When  the  quasimode  is  negative  energy  in  the 
beam  frame,  the  beam  transfers  energy  to  it.  When  the  electro¬ 
magnetic  wave  is  negative  energy,  energy  flows  from  EIC  to  EM 
waves  preferentially. 
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EIC  wave 


X-mode 


beam 

Fig.  9.  Another  possible  version  of  the  mechanism  to  produce  kilometric 
radiation.  The  beam  naturally  radiates  an  R-mode.  This  mode 
provides  a  conduit  for  EIC  waves  to  feed  energy  into  the 
negative  energy  X-mode,  which  then  grows. 
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